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We investigate the basic theoretical issues in the quantum entanglement of particle pairs created 
from the vacuum in a time-dependent background field or spacetime. Similar to entropy generation 
from these processes which depends on the choice of physical variables and how certain information 
is coarse grained, entanglement dynamics hinges on the choice of measurable quantities and how the 
two parties are selected as well as the background dynamics of the field or spacetime. We discuss 
the conditions of separability of quantum states in particle creation processes and point out the 
differences in how the von Neumann entropy is used as a measure of entropy generation versus for 
entanglement dynamics. We show by an explicit construction that adoption of a different set of 
5^ ' physical variables yields a different entanglement entropy. As an application of these theoretical 

Oh' considerations we show how the particle number and the quantum phase enter the entanglement 

, dynamics in cosmological particle production. 

O ' PACS numbers: 03.65.Ud, 03.65.Yz, 04.62.+V 

CO 



o 
cr. 



CN 
> 



X 



I. INTRODUCTION 



Entanglement is said to be the uniquely distinguishing feature of "quantumness" [l|. Despite increasing interest 
in the recent decade and advances in quantum information sciences, we are still far from fully understanding the 
nature and dynamics of quantum entanglement - how it is characterized and how it evolves in time - in quantum 
^vq ' open systems, i.e., those interacting with their environments. We believe at this stage of development it is important 
£SJ \ to analyze these basic issues in great detail in simple enough systems, preferably with exact solutions, so that we can 
0^ ■ understand in depth its behavior to gain much needed valuable insight 

Usually entanglement is discussed in the framework of nonrelativistic quantum mechanics, but for a quantity to have 
_^ . physical meaning one needs to know how it transforms in different reference frames, e.g., for two observers moving at 
relative constant speed how is the quantity one observer reports as entanglement between two parties in its system 
related to that reported by the other observer in its system? To answer this rather rudimentary question one needs to 
work with Lorentz transformation of entanglement in the framework of relativistic quantum mechanics (see, e.g., Q 
for a review). This is the starting point of relativistic quantum information (RQI). When a quantum field is involved, 
one needs to upgrade the treatment to that of relativistic quantum field theory. This is the platform we have adopted 
for our research program on RQI. When quantum informational issues arise pertaining to black hole information loss 
or early universe quantum processes we need to consider them in the extended framework of quantum field theory in 
curved spacetime 

The simplest process which distinguishes a quantum field theoretical process from a quantum mechanical one 
is particle creation. The investigation of quantum entanglement in the particle creation process - how to define 
entanglement, between what parties, and how it evolves in time - is the first order of business toward establishing a 
RQI theory for quantum field processes. This is the goal of this paper, focusing on particle creation in strong and 
dynamical background fields, such as in the Schwinger effect Q and in background spacetimes, as in cosmological 
particle creation [l(J [H|. ^ n addition to its theoretical value these results are expected to be useful for quantum 
information experiments and for probes into the very early universe from next generation cosmological observations. 

The statistical mechanical properties of particle creation such as entropy generation have been a subject of both 
theoretical and cosmological interest for quite some time. Since the mid-1980s there are inquiries on finding a viable 
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measure of entropy for particle creation processes in a free quantum field [12J and from interacting quantum fields 
fl3| . (The most detailed discussion of entropy generation from free field particle creation along the theme of [l2[ can 
be found in [lij]. For a summary of recent developments see Chapter 9 of 15].) The free field case is conceptually 
more intriguing: On the one hand, from a pure quantum field theoretical perspective, one would say categorically 
there could be no entropy generation because the particle pair is originated from the vacuum which is a pure state. 
On the other hand, from thermodynamical considerations one sees clearly that entropy is generated in the amount 
proportional to the number of created particles. This puzzle is what started the inquiries mentioned above which 
led to the understanding that one could quantify entropy generation by the number of particles created only because 
one chooses to measure this process in the Fock (number) representation at the sacrifice of the phase information, or 
effectively imposing a random phase approximation. One also knows that it is only for the spontaneous production of 
bosons that particle number increases monotonically, a necessary condition to associate it with an entropy function 
obeying the second law. For fermions and for stimulated creation processes particle number can decrease which 
invalidates this definition. What this earlier investigation taught us is that the measured entropy associated with the 
particle creation process depends on quite a few factors: the choice of physical variables (such as using the number 
basis), the coarse graining one introduces in its measurement (such as the quantum phase) or recognizing that some 
information is inaccessible to the observer in the process. 

We mention prior studies on entropy generation in particle creation processes because the experience we gained 
(concepts and methods) and the lessons we learned (see above) prove to be useful for our investigations into the 
quantum entanglement of particle creation because they share certain qualitative similarities. Specifically, quantum 
entanglement depends crucially on the choice of physical variables and the way the two parties whose entanglement 
is the object of interest are defined. Quantum entanglement can be measured in many ways [Tr| . For the bipartite 
system the von Neumann (vN) entropy is quite commonly used. We will show that the way the vN entropy is used 
and what results it yields for these two processes, one pertaining to entanglement dynamics, the other for entropy 
considerations in nonequilibrium statistical mechanics are quite different, both in terms of whether and how coarse 
graining is introduced, and for what reasons. A particle pair with some 3-momenta (k, — k) created from the vacuum 
is habitually regarded as perfectly entangled. We show by explicit construction that adoption of a different set of 
physical variables with the same 3-momenta (k, — k) makes a difference in the entanglement dynamics. 

Before we go into detail, let us begin with a description of how information is chosen or coarse-grained in a closed 
versus an open quantum system. (For a treatment of entropy generation in squeezed open quantum systems with 
application to inflationary cosmology, see [1 71 ] and related work cited therein.) Any quantum state describing a closed 
(entirely isolated) system is always pure, and the vN entropy is identically zero. In realistic settings a system is 
rarely closed as it interacts, no matter how weakly, with its environments. The back action of its environment with 
some coarse graining (in which process noise is engendered) induces dissipation in the dynamics of the system, now 
rendered open, and mixed states are allowed in the open quantum system. This is one way entropy is generated in 
the system. Another way this could happen is either by necessity, that some degree of imprecision always exists in 
realistic measurements (or uncertainty in quantum state tomography), or by choice, that only some physical variables 
of direct interests are measured (e.g., number rather than phase in particle creation) and others ignored, averaged, 
or "integrated out" (e.g., the fast-oscillating elements of the density matrix). These procedures either by necessity or 
by choice would render a pure state mixed in appearance and entropy generation ensues. 

The first type of entropy generation is well illustrated in an open quantum system treatment of cosmological particle 
creation such as is shown in [17] , and more powerfully in the particle creation and backreaction problems, where the 
gravitational sector is viewed as the (open) system and t he q uantum field as its environment [l|| , whose backreaction 
causes the dissipation of anisotropy or inhomogeneities |l9l . |20| of the early universe. One can obtain the entropy 
gen erated in the particle creation process in terms of a vacuum viscosity function, even define gravitational entropy 
[21| of spacetime dynamics associated with these processes. The second type is illustrated clearly in free field particle 
creation processes [HI, [lij , where the choice of a Fock representation is what enables one to relate the amount of 
entropy generated to the particle numbers created. Considering particle production in a uniform electric field Kluger, 
Mottola, and Eisenberg (KME) 14] derived a quantum Vlasov equation describing the evolution of the (adiabatic) 
particle number, and an equation for the evolution of the quantum phase. They then argued that in the case they 
considered since the quantum phase as well as the off-diagonal elements of the density matrix (the coherence) oscillate 
rapidly in time, these fast variables are expected to be averaged out in observations and the density matrix would 
effectively look like a mixed state to such an observer. 

Oftentimes these two sources of entropy generation are intermixed, depending on how a problem is formulated 
and treated. This happens in the interacting field particle creation calculations of [l3j |. A more recent example is 
in [22l [23[ , where Campo and Parentani studied the self-consistently truncated "Gaussian and homogeneous density 
matrix" (GHDM) for an interacting field. They calculated the vN entropy of the two-mode sector with opposite wave 
vectors (k, — k) of the GHDM and showed that the vN entropy is the only intrinsic property of the field state during 
inflation, while the entanglement between particles with k and their — k partners depends on the choice of canonical 



3 



variables in the same mode pair </>k</>-k- Their vN entropy of the two-mode sector plays a double role. On the one 
hand, it is a consequence of truncation or coarse graining in obtaining the GHDM, effectively by averaging out fast 
variables. On the other hand, it is a measure of the entanglement between that mode pair </>k0-k and its environment 
consisting of all other modes (rather than the entanglement of k and — k particles in the same two- mode sector), by 
noting that the GHDM is factorized and the two- mode sector with momenta (k, — k) is by itself the reduced density 
matrix (at least approximately) with other degrees of freedom traced out. 

To get rid of such intermixing, we are focusing on the simplest cases with free quantum fields in a dynamical 
background. In Sees. II and III we describe the particle creation process of real and complex scalar fields, respectively, 
in the Schrodinger representation. In Sec. II B and Sec. IV we investigate the behavior of particle numbers and 
quantum phase exploring the theoretical issues for the quantum entanglement of particle creation in a time-dependent 
background. Section II C contains the main results of entanglement dynamics using Wigner functions. After the 
theoretical issues are explored and analysis performed we study such processes in the early universe in Sec. V. We 
find that for the vacuum state of a free real scalar field in an expanding universe, once the physical variables are 
correctly chosen, it is possible to partition the degrees of freedom of a (k, — k) mode pair into k and — k particles, 
and the degree of entanglement between them can be calculated accordingly. Based on these results we are able to 
look into how the particle number and the quantum phase enter the entanglement dynamics in cosmological particle 
production. We conclude in Sec. VI with a brief summary of the key results reflecting the main themes stated here. 

II. REAL SCALAR FIELD WITH TIME- VARYING MASS 

Without loss of generality, let us work with a free scalar field in Minkowski spacetime with a mass parameter M(rf) 
which is time dependent, representing all time-dependent parameters entering into the system (including that from 
cosmological spacetimes) 



S 



(1) 



where we denote time x° by r\ anticipating cosmological applications later. In Fourier-transformed representation 
4>(x) = (27r)~ 3 J c? 3 /ce lk x (/) k , the above action becomes 



s= dn 



d 3 k 



~<3,,<M^-k " T^kM^k^-k 



(2) 



where ^£(77) — k 2 + M 2 (rj) and </>_ k = <f> k . The momentum conjugate to <p k is Ilk = SS/S(dr,(pk) — c^-kj and the 
Hamiltonian can be derived straightforwardly by a Legendre transform of S. 

In Sec. IV Al we will see that the theory with a real scalar field in a Friedmann- Robertson- Walker (FRW) spacetime 
has the above form in conformal time rj. 

A. Quantization in Schrodinger representation 

Canonical quantization of this theory is achieved by imposing the following equal-time commutation relations, 

= ^) 3 <5 3 (k- P ), 
[faWMv)] = [n k (»?),n p (r ? )] = o. (3) 

In the Schrodinger representation (24|, ^> k are viewed as c- number functions and the conjugate momentum operators 
are given by 



nk = (27r) 3 ^, 



(4) 



to satisfy the above commutation relations. The Hamiltonian operator then reads 



H 



(X) 3 



i 00k 0<P-k ^ 



(5) 
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In this representation all physical information of the system is included in the wave functional ^[4>\.,rj\ satisfying the 
Schrodinger equation 

ihdrfS = H^. (6) 
The ground state with minimum ( H ) is given by the Gaussian state 



-' rw<,)/ *->4/(0/|?^^ w 

where s is the normalization factor and Eq is the vacuum energy [361 ]. Substituting the above ansatz into ([5]), one 
finds that 

.9 k ' p W = ^4rT( 27r ) 3,53 ( k + p) ( 8 ) 



so that = Ilk^ok' w here 



* 0k = * k e J ^W/ exp ____^_ k; (9) 
with the normalization factor Sk and the "ground state energy" for the k or — k modes, 

£ k ^[lx k | 2 + fi k |x k | 2 ], (10) 

such that E (ri) = Ek^o- Here J2 k = J dkd(0), a prime denotes taking the derivative with respect to ry, and the 
mode function amplitudes x k satisfy the equation 

Xk + ^(*7)Xk=0, (11) 

so that Xk = X-k- This as we recognize is the classical field equation for the scalar field, and \k 1S the solution for 
this field equation subject to the normalization condition 

XkXk ~ XkXk = * (12) 

at every moment, which requires Xk to be complex functions. Note that Xk does not have to be classically physical 
solutions, because boundary conditions for Xk wm not be fixed by classical considerations. In fact, different Xk could 
correspond to different quantum states. Compared with the ground state for real scalar fields with constant mass 
in Minkowski space (26|, one sees that Xk should be taken as \/l/2f2fce _^0 ' c ' , there, which is the positive frequency 
component of the kth mode of the free scalar field. 

The wave functional in the form Q with Xk the solutions of the classical field equation is almost the same as an 
infinite product of the wave functions for squeezed Gaussian states of parametric oscillators with frequency £lk{r]) 
[25j . except that the variables in the exponents are mode pairs (f>k<fi-k for each k rather than the squared "positions" 
of the oscillators. 

The nonvanishing two-point correlators in the vacuum state (J7J are 

(0 k ,0 p ) = M2^) 3 <s(k + P )| Xk | 2 , (n k ,n p > = M2^) 3 <5(k + p)|x k | 2 , (13) 

and ( cj) k , Hp ) = ^drj ( <p kl (f) p ), where ( A,B ) = | J T>(f>^*(AB + BA)^ . Then one can write down the covariance 
matrix V with elements Vy = (lZi,lZj ) in 7£j = (</> k , II k , 0_ k , n_ k ) representation, whose partial transposition 
V PT = V|n_i,->— n_i, gives the quantity we used in an earlier paper [3, H3] 

£ = det^f T + i[^,% 



i[M2^ 3 (o)] 4 [(a r) | Xk l 2 ) ; 



(14) 



Since this is always positive, the field modes </> k and 0_ k are not only separable from other degrees of freedom in the 
vacuum state (i.e., the von Neumann entropy of the density matrix l I / o k , I , Q k is zero, or in other words, each ^ok is 
always pure and factorizable from ^0)1 but also "separable with each other" at all times. 
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To see this more clearly, we write the complex </> k in terms of two real field variables as k = </> k + i<p^ instead, 
with </>^ k = k and = — c/f> k (since = 0_k)- Then one can easily see that ([§]) can be factorized into a product 
of k state and ^ k state. So if the observables are (</> k , k ) fields, they will always be separable with each other and 
no entanglement measured in terms of these variables will be generated for a free scalar field with time- varying mass, 
such as the field in an expanding universe. 

Nevertheless, quantum entanglement depends on partition as well as the choice of physical variables or measurables. 
With reference to quantum entanglement, say, in cosmology, foremost one needs to specify which physical observables 
are being measured there. Obviously k and k are not the correct variables to describe quantum entanglement in 
cosmological particle creation: Recall that modern cosmological experiments measure the temperature fluctuations 
6T/T of the cosmological microwave background radiation rather than the amplitudes of quantum fields. The former 
can be related to the energy density perturbation Sp/p with the energy density p = (H) ~ J^k { Nk ) /V, where 
H and iVk are the Hamiltonian density and the number operators for <fi k , respectively. The number operator 
here consists of creation and annihilation operators defining the in/out vacuum at the initial/final moment (or in 
the adiabatic vacuum). We therefore intend to look at quantum entanglement in terms of these operators. This also 
suggests that quantum entanglement generation, like cosmological particle creation, manifests only in those physical 
variables which facilitate a well-defined in/out or adiabatic vacuum. 



B. particle numbers 



Continuing our exposition using quantum field theory (QFT) in Minkowski space, we define the annihilation and 

k( 



creation operators & k (r?) and b^_ k (rf) by 



b ^ = /*,o ^ 3 ^ (Xk*(^k-Xk(??)n-k), (is) 

so that 

k = #^M0) (xu(v)b k (v) + xi(v)bU(v)) , (17) 

n_ k = VM2^(0) (xL(v)b k (v) + Xw(v)bl k (v)) , (18) 

are independent of time. The above definition of operators fr k (?7) and b_ k (r)) has the following properties: First, 

they become the conventional ones in QFT in Minkowski space where x k = y/l/2ilks~ ink n ■ Second, the equal-time 
commutation relations Eq. (J3J) are equivalent to 

[6 k ,6t] = j3 (k _ p)/( j3 (0)< (19) 

Third, the ground state is the state with the property fok^o — for all k. Fourth, the "excited states" are 
analogous to those in simple harmonic oscillators, which are generated by applying these creation and annihilation 
operators to the ground states, e.g., 

*> (k) = b --*°= vmFm k*°' (20) 
..o, - k) - ^u*. - (s^l^ - 1) ••■ (21) 

are also solutions of the Schrodinger equation. Moreover, all excited states generated in this way, together with the 
ground state, form a complete set of quantum states for the scalar field. 

Note that a proper normalization for 6 k and 6_ k has been chosen to make the above excited states satisfy the same 
normalization conditions for vPq. Written in terms of 6 k and fr_ k , the Hamitonian operator reads 



H = E { £ o (bA + &&) + ^ [(X k 2 + «kXk 2 ) M-k + H.c] } 



(22) 
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where "H.c." stands for "Hermitian conjugate" and £q is the kth component of the ground state energy given in (|10l) . 
So J T>(f) '^Q^H'i 1 (r]) = E (r]) is indeed the vacuum energy. 

The "number operator" of mode k at the moment rj is defined as 

MV) = bt(v)b k (v)- (23) 

Suppose at the initial moment 770 the field is in the vacuum state when the particle-number counter is constructed 
according to the operator N^(r}o). Then, using the same particle-number counter, one finds that the number of the 
particle created at time 77 by the background spacetime is 

(O v \Nu(Vo) K) = J ■D4>%(r ] )N k (Vo)%(v) = \Xk(vo)Xk(v) ~ Xk(vo)Xk(v)f ■ (24) 
Since both XkiVo) an d Xki 7 ]) are solutions of and each form a complete set at different times, one may write 

XkM = a k ^o,r])x k {v) + Pk(VQ,v)xUv), (25) 
or equivalently, for |a k | 2 — \(3 k \ 2 = 1, one has 

Xkiv) = (xlXkiVo) - AcXkM- (26) 

Then the Bogoliubov coefficients read 

"k = i[Xk(m)Xk(v) - XkMx'kiv)} +i((v)xl{ri)x k (Vo), (27) 
Pk = i IXkMx'kiv) - XkMXkiv)} -K(v)Xk(v)Xk(Vo), (28) 

with overdots denoting d/dr]\ v= , l0 and (e R, whose values will be further fixed by specifying x'k or other physical 
conditions. From ([24| it can be seen that we have made the physical choice £ = 0, which implies (0^| Nk{vo) |0^) = 

Note we did not get into the details about how to define an adiabatic particle-number state (which can be found 
in e.g., [1, HH) but simply assume that at rjo the concept of particle is well defined. Here the Bogolubov coefficients 
with £ = are represented explicitly by the Klein-Gordon inner product. 



C. Wigner function and entanglement entropy 

Since \I/o is a pure state, if one divides the degrees of freedom in this model into two parties, then the entanglement 
between them can be well measured by the von Neumann entropy of the reduced density matrix of one of the two 
parties. Nevertheless, it is not obvious whether the particles with k and — k are separable as indicated by the positive 
E at the end of Sec lTTAl 

A clearer separability can be seen in the Wigner function of ^>q: 



** l,r/TJ3 fill ^ L \ / \ ~ 



ke(R3-{o})/z 2 

2 f d 3 k 

1 pvn / 

hj (2tt) 



»\'^V-T I ^{iXkWl'^-k + IXkWl'^-k-a^lXkWI 2 )^^}, (29) 



where we have taken <^ k = + i<f) k (with </>^ k = k , = — (j> k ) and ir±k = (p k ± ip k )/2, which are c-numbers 
here rather than operators. Performing a canonical transformation similar to flTTl) and (TT5|) as 

^ = v/ft(27r)3«S3(0) (xk(v)Bk(v) + xi(v)B\(v)) , (30) 

n k = #^M0) (xi(v)Bk(v) + Xk*(v)B*-k(v)) > (31) 

which gives dcf> k dir k — dB k dB k h(2Tr) 3 5 3 (0) , one ends up with the Wigner function in a Fock representation, 

p{r,) = |5| 2 exp-2^ J B k (r,)B k (ry). (32) 
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Here one can easily see that the quantum states of (B k , B k ) and (B_ k , B*_ k ) for each specific k are separable. 

However, in terms of an alternative set of measurables, k and — k particles could be entangled as we will see in the 
following. From (|24|) it seems that the physical particle counter should be counting N^tjo) defined by operators 6k 
at the initial moment r)o. Let B k = B k (r]o) in (I3TJ1) and (f5T|) at r] = rjo, so that 

B k ( v ) =« k i? k + /3 kJ B* k , (33) 
from (f2l)]) . Then in terms of B k with nk = |Ac| 2 , Ck = «k/3k, the Wigner function (|32p becomes 

p( V ) = |S| 2 exp -2 £ [(2n k + 1) B k B* k + c k B k B_ k + c* k B* k B*_ k ] , (34) 

k 

which gives the "symmetric" particle number 

( B* k B k ) = n k + ~ = (0„| fet^o), 6 k M |0„) . (35) 

This justifies B k as the correct variables corresponding to physical measurements. Suppose we divide the particles 
into two groups with fc 3 > and fc 3 < respectively (that the particles with fc 3 = are properly divided into these 
two groups is understood). One may write 



p( V ) = |S| 2 exp -Y, + [(4» k + 2) (B k 5k + S-k^* k ) + 4c k B k B_ u + 4c k £ k i?* k ] • (36) 

k 



where ^2 k + an< ^ l a t er Ilk + denote summing and multiplying over k with fc 1 , fc 2 6 R and fc 3 > 0, respectively. 
Integrating out B_k and B*_ k in p(r]), one obtains the reduced Wigner function 



P R (V) = I n +dB -k^* k p(r,) 



| s -| 2 exp-]T + |^[(2n k + l) 2 -4|ck 

k 



2n k 

k 

which is a mixed state once |/?k| 7^ because the purity 

7> = 2tt I dB k di? k (/) 2 = I] + (38) 

k 

is less than 1 if any n k = |/?k| 2 > 0. This means that the particles measured by "physical particle counter" iVk(»7o) 
with k and — k are entangled. Comparing the above expression for the purity with Eqs. (30) and (31) of Ref.[f|, one 
sees that the von Neumann entropy is 



S = 



^ [(ra k + l)ln(ra k + 1) -?i k lnn k ] 

k 

= ^[( nk + 1 ) ln ( nk + 1 )~ nklnnk ]- ( 39 ) 
k 

When n k >• 1 for all k, S ps l nri k- 

vN entropy as an entanglement entropy measures the quantum nonlocal correlations at some moment between the 
system and some specified party which could be the rest of the world. The latter is traced out at the moment the 
entanglement entropy is evaluated. Note that "nonlocal" here does not imply "nonlocal in space": the particles in 
this paper are actually something similar to plane waves, which are local in momentum space but nonlocal in position 
space. They are not spacelike separated. Rather, they can occupy the same position space at the same time. 
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III. COMPLEX SCALAR FIELD WITH TIME- VARYING MASS 

One can easily generate the above formulations to a complex field with time- varying mass: 

1„ 1 



S = J d 4 x --d^d"®* --M 2 ( V )\$\ 2 , (40) 
or in a Fourier-transformed representation, 

~ar„\\A 12 ( 41 ) 



(2tt) s 



The only difference is that now is independent of 0_ k , while II k = d v (j) k . The above action describes scalar QED 
of a quantum charged scalar field in a classical uniform electric field [l4T |. 

The ground state wave function is still in the form of ((9]), except that the 4>_ k should be replaced by But now 
the particles and antiparticles can be distinguished by their different charges in addition to their momenta. So instead 
of (JT7J) and HH]), one writes 

^ = VW(0) (x k (*?K(»?) + xtfa)&l k fa)) , (42) 
n k = V^(2^) 3 ^ 3 (0) (xLWokW + XkW&tkfa)) > (43) 

with 

[a k ,at] = [6 k ,6t] = J 3 (k-p)/^(0). (44) 
This will make the Wigner function of the vacuum state a product of two copies of ([32]) : 

pfa) = |S| 2 exp ~2 J] (l k (r?)l k (??) + S k (r?)B k (r?)) , (45) 

k 

where (^4. k ,-B k ) corresponding to (a k ,& k ) are the counterpart of B k of the real scalar field. Thus as before one can 
express them in terms of the ones defined at the initial moment as 

A k ( V )=a k A k + p k B*_ k , B_ k (r,)=0 k A k + a* k B*_ k , (46) 

according to (|25|) . These imply 

P ( V ) = \S\ 2 exp -2^ [(2n k + 1) (A k A* k + B k B* k ) + c k (A k B_ k + A_ k B k ) + c k (A k B*_ k + A*_ k B k )] . (47) 

k 

Although a mixing between the particles with k and the antiparticles with — k is generated, the vacuum state ^0(77) 
remains a pure state. Only after one integrates out the antiparticles (particles) associated with B k (A k ) will the 
reduced Wigner function for particles (antiparticles), p A (p B ), become 

^) = | s -| 2 exp-]T|^, (48) 

k 

with C = A,B. Now p c is a mixed state. The purity and the von Neumann entropy are those for the real scalar field 
<[3H|) and (f3T))) with Y[ k and J^u replaced by the usual Y\ k and J^ k . So the values of the von Neumann entropy of a 
complex scalar field between particles and antiparticles are twice the value for a real scalar field, indicated by (|39|) . 

IV. PHASE INFORMATION 

Observe that the description from (|32p to (|3"4")l is that of squeezing in a two-mode squeezed state, well known from 
a squeezed-state description of particle creation (see, e.g., [13, HH). Writing B k in terms of quadrature amplitudes, 
namely, B k = (Qk + iP k )/v2 with Q k and P k real, then <f32]) looks like a direct product of the Wigner functions for 
the ground states of Q k for all k. Now, since |a k | 2 — IAc| 2 = 1, one is allowed to parametrize a k and (3 k as 

a k = e mk coshr k , f3 k = -e i0k sinhr k . (49) 
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Let B k = (Q k + iP k )/\/2. Then one has 

Qk = cosh r k ^Qk cos Ok + -Pk sin a\^j + sinh r k ^Q-k cos k — P_ k sin 9\^j , (50) 

P k = cosh rk ^—Qk sin CTk + Ac cos f7k^ " sinh rk ^Q-k sin ^k + P-k cos 0k^) , (51) 

from the relation (j3"3")) . We see this involves two steps: 1) (Q±k,-P±k) are first rotated locally in angles er k and — # k 
on the +k and — k modes, respectively, and then 2) squeezed globally in squeeze parameter rk- Local operations do 
not affect the entanglement measure: The phases cr k and k are not present in the reduced density matrix (RDM) 
transformed from (|37l) , neither is any function of the RDM such as the entanglement entropy (1391) . Here only squeezing 
which is a global operation is relevant to quantum entanglement. 

We now turn to the question of how to obtain the phase information and whether/how it enter into entanglement 
dynamics considerations. 



A. Quantum phase 

Mathematically the phase # k of /3k in the parametrization (l4"9l can be obtained by evaluating 

#k 



and Ok in «k can be obtained in a similar way. But physically only the phase sum Ok + #k could be measured. 
Transforming Q34[) to the density matrix in quadrature amplitudes Q± k representation, one obtains 



p[- ■ ■ , Qk> Q-k> ' 



, Qk: Q-k, •••]=*[..., Q k , g_ k , • • .]*•[. ■ • , Q' k , Q'_ k , • • •], 



(52) 



(53) 



where * = Ilk^k and 



e -2f v dfj£*(fj)/h _ x 

*k = — exp— {(l + 2n k + c k 2 -c k 2 ) (g^+Q^) +4 [c k (n k + l) -c k n k ]Q k g_ k } 



cosh rk 



1 + sinh 2 r k (l - e 2i ^ k+e ^) 
1 + sinh 2 r k (l - e- 2l (^+ e ^) 



1/4 



£ (j^+V tanh r k ) " $„ (Q k ) $„ (Q_ k ) , (54) 



n=0 



with (? = (! + 2n k ) — (c k + c k ) and the number eigenstates in Q representation, 



So the density matrix (|53[) can be expressed as 

P = Il + E (Qk) (Q-k) *™ (Qk) *m (Q-k) 

k n,m 

with 



k = tann r u i(n _ m)( ^ + g y) 

f nn.mrn i 2 

COSH 7Y 



(55) 



(56) 



(57) 



Equation (|54l) shows that k and — k particles associated with Q k are always created in pairs, because the outcome of 
the measurement on numbers of k and — k particles separately will always be the same. From (|56[) and according to 
[29l | , one can write down the probability distribution of the quantum phase-sum 9 + , 



P(9+) = {2-7T [cosh 2 r k + sinh 2 r k — 2 coshr k sinhr k cos (9 + — er k — k )l } 1 , 



(58) 



which peaks at 9 + — ((T k + # k ) = 2mr, n € Z, while the probability distribution of the quantum phase-difference 9-, 
P(9-) = l/27r, is constant through < < 27r so 0_ of the quantum state (|5~4"|) is totally uncertain. This shows that 
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[(<7 k + # k ) mod 2ir] is the quantum phase complementary to the particle number, while Ok and 6*k cannot be observed 
separately. 

The matrix elements (|57j) are proportional to c k _m [37(. If the phase changes so fast that c k could not be measured 
precisely by any apparatus, then the information for quantum state tomography will never be complete. Most likely 
in this case the off-diagonal elements with m ^ n would be averaged out, then the quantum state of the field may 
appear like a classical state. This "fake decohcrcncc" due to the technical limitation of measurement is different from 
environment-induced decoherence such as from intermode couplings. 

If the resolution of the apparatus gets higher, more phase information could then be observed. However, it is easy 
to verify that once p^ n mm oc c k _m for all nonzero c k € C, the vN entropy of p will be zero. So vN entropy cannot 
be generated by just replacing all the original c k by some outcomes of measurement with smaller absolute values. 
One simple way to produce vN entropy is to perform a truncation such as p\ n mm = for all \n — m\ > N with 
some positive integer N, meaning that the off-diagonal elements oscillating quicker than exp±iiV(cr k + 6* k ) are not 
resolvable and being averaged out. For this truncated effective density matrix, the purity will be 



V, 



2n k 



2 



N 

I Tit. 



^ \ 1 

m— 1 



(59) 



For all nk > 0, the larger N is, the closer V e s is to unity, so the effective density matrix is purer, and the vN entropy 
of it is closer to zero. 

The RDM obtained by tracing out the Q_ k and Q'_ k components in (151))) reads 



r 

k 



P R = Il + E Pnn,nnK (Ok) K (Ok) ■ (60) 



One can immediately see that the vN entropy of (|60l) . which is the entanglement entropy between the particles with 
k and — k, has exactly the same value as the vN entropy of the effective density matrix of the vacuum with all 
off-diagonal elements averaged out, namely, 

^cff = II E Pnn^n (Ok) (O-k) K (Ok) K (O-k) ■ (61) 



Thus one could say that the tracing-out process in obtaining the RDM (|60[) represents the "strongest" coarse graining, 
though the coincidence of the entropy values here does not occur for general quantum states. 

One should be careful that simply ignoring fast-oscillating elements could create one more problem if our quantum 
state tomography is designed to reconstruct the Wigner function. The amplitude of Ck plays an important role in 
obtaining the correct entanglement entropy from the Wigner function. If one finds that Ck in (1341) appears to be zero, 
then the factor in the exponent of the reduced Wigner function ([57]) will be 2(2nk + 1) rather than 2/(2?ik + 1), so 
that the effective reduced Wigner function cannot be transformed back to the correct effective RDM. This effective 
reduced Wigner function yields an incorrect vN entropy or purity for the k particles, though here the vN entropy is 
no longer a well-defined entanglement entropy since the density matrix of the (k, — k) mode pairs (not the RDM of 
the k particles) effectively constitutes a mixed state. 

Since the quantum phase is conjugate to the particle number of a squeezed state, one may expect that one could 
obtain the phase information from time derivatives of nk or equivalently, from evolution of the entanglement entropy 
in time. Indeed, by noting that d m Xk(v) = an( i lrom p6])- one has Xhiv) — ^(Vo, w)Xk{Vo) ~ Pkivo, v)xt.( r lo)- This 
implies that 

< = -2 Im {c k [x k 2 (%) + ^)Xk 2 (?7o)] } , (62) 

where c k = a k /? k = — e J ( <Tk + 6 ' k ' coshrksinhrk provides information of the phase <7k + #k- Unfortunately, in the right 
hand side of (|62[) . c k is always multiplied by a term in the square bracket, which often cancels the oscillation of c k so 
that one cannot read off the quantum phase from the behavior of n k . 

Even c k per se are not always fast oscillating, though. Below we will give an example when the off-diagonal elements 
of p\ n mm associated with the particle-number operators defined at the initial moment do not oscillate, namely, when 
the Universe undergoes inflationary expansion. But before we proceed, we want to make one more remark on an 
alternative oscillating "phase." 
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i(2n k 



B. Quantum Vlasov equation 

The r\ time derivative of c k reads 

4 = 2ic k [|x' k (%)| 2 + r! 2 (r,)|x k (r/o)| 2 " 
Similar to iTi} , one can express Ck in terms of nk by solving (|63|) then insert it into (|62j) to get 

n k (r?) = 2Re { [ Xk 2 (vo) + ^l(v)x k 2 (vo)] x 

' dr? (2n k (ij) + 1) [x k 2 (%) + f^M'M] e^ 9 ^)" *^ 



where 



e k (r?) 



IXk(%)| 2 + ^(r)|x k (%)r 



(63) 



(64) 



(65) 



which demonstratively illustrates that the evolution of n k is in general non-Markovian. Only in the case that the phase 
| © k (77) — @ k (77) I grows rapidly in 77 — 77, would the fj integration be effective only around fj « 77, and the right-hand 
side of (|64]l becomes local in time. 

Anyway, solving the quantum Vlasov equation (|64[) . or (|62[) and (|63p . is equivalent to solving % k from (|11|) and 
then calculate n k and c k , which is much simpler. We put (|62[) and (I63|) here simply to show the relation between the 
particle number and the phases. Often it is not economic to solve them directly. 

We should emphasize that O k in ([65|l and the counterpart in [lj] , which is the phase of the adiabatic mode function, 
are different from the quantum phase cr k + # k in general. 



V. ENTANGLEMENT IN COSMOLOGICAL PARTICLE CREATION 



A. A real scalar field in the FRW universe 



A real scalar field $ with mass m minimally coupled to a curved spacetime with metric g^ is described by the 
action, 



S = J d^x^—g 



2 2 



(66) 



We are working with a test-field condition where the gravitational field g^ is a given background, in this case, the 
FRW spacetime, with line element 



ds 2 = a (77) 5 



-drj 2 



dr 2 



1 — nr 



■ .2 



r 2 d0 2 + r 2 sin 2 



a(rj) 2 \-dr/ 2 + hijdx % dx jl \ 



(67) 



where k = 1,0,-1 corresponds to closed, flat and open universes respectively. In terms of the conformal time 77 and 
the conformal scalar field defined as 



X(x) = a(ri)$(x), 



one can rewrite the field action S as 



-X' 2 - -diXd'X + -(—- m 2 a 2 ) X 2 
2 2 2 V a 



S = J dr]d 3 xVh 

plus a surface term — / d 3 xX 2 a' /2a. The field equation then reads 



X" + m 2 a 2 



V A = 0. 



(68) 



(69) 



(70) 
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To better handle the spatial derivatives in action S, we perform a transformation 

X(:r)=^ k (x)</. k (77), (71) 

k 

where and 3^c(x) for k = —1, 0, and 1 can be found in [8]. For n — (spatially flat), 5Zk' = J c? 3 ^/ (27r) 3 and 3^k 
is simply e lk x , and the action written in the Fourier k space becomes (J2]) with time- varying squared frequencies 

!l^)Ei 2 |mV-^. (72) 
So the formulation in Sec. [II] can be directly applied. 



B. Entanglement of particle creation in a de Sitter spacetime 



In the spatially fiat FRW coordinatization of the de Sitter space, the scale factor is a 



-(Hri)- 



with 



Hubble constant H and cosmic time t, such that rj runs from — oo to as t goes from — oo to oo (see, for example, 
Sec. 5.4 in Ref.[8j]). The squared time- varying natural frequency in (|72p now reads 



nl = k 2 



m 



1 



(73) 



If m 2 < 2H 2 , Sl k will become negative at late times (rj 
to taking the value [30| 



0) for all finite k. The Bunch-Davies vacuum corresponds 



x k = yfi4 2) (M, 



(74) 



where is the Hankel function and v = [(9/4) - (m 2 / H 2 )] 1 / 2 . 

The phase in the quantum Vlasov equation (I64[) of such a scalar field in this de Sitter spacetime 



^(v)-Qu(v) = (v-v)<\xM\ 



k 



H 2 



- 2 



1 

VV 



ixkMr 



(75) 



varies quite rapidly in 77 — fj time, so the integration in (|64p is more pronounced around fj ~ rj and the behavior 
of n k in an inflationary universe can be treated in a Markovian approximation. However, this does not imply that 
the quantumness of the field is lost in this epoch, as some authors working on the decoherence of quantum fields in 
inflationary cosmology are drawn to making such a claim. 

Substituting (|74[) into (|52p and the counterpart for er k , we find that cr k + (9 k varies quite slowly except when 
rj 2 \[{m 2 / H 2 ) — 2]/k 2 \, which is around the moment that f2 k is crossing zero in the cases with m 2 < 2H 2 . Far 
from this moment in cosmic time t, c k _m almost does not oscillate and could be identified clearly for very large 
\n — m\. Therefore perhaps contrary to common belief, in almost the whole inflation epoch, the off-diagonal elements 
of the density matrix corresponding to the quantum interference between the k and — k particles associated with £> k 
corresponding to the in- vacuum manifest and the coarse-grained effective density matrix is extremely pure j38| . 

Alternatively, in the adiabatic number basis, the two concepts of phase in Sec. IIVI can coincide. For example, the 
(first-order) adiabatic mode function used by KME in [l4[ reads 



1 



2ftk(f?) 



exp 



n 

-1 I d?7f2 k (?7). 

'no 



Numerically we find that at early times when f2 k is real and not very small, the adiabatic particle number [l l 

Wk = |xk (Xk + *^kX k )| 2 



(76) 



(77) 



is indeed much less than n k , and the counterpart of c k for the adiabatic mode function does oscillate, in exactly 
the same way as the oscillation in the phase in their quantum Vlasov equation (which implies that it is still 
impossible to determine the quantum phase by observing the evolution of their A/"k)- This justifies the argument in 
[14J : The adiabatic particles look more classical since the off-diagonal elements of the density matrix of the vacuum 
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(the counterpart of (|57[) with r k and <7 k + k obtained in the adiabatic number basis) oscillate too fast to be resolved. 
Then the vN entropy S e g of the effective density matrix p cS of the vacuum with all off-diagonal elements in the 
adiabatic number representation vanishing is (|39j) with n k replaced by Afk [cf. Eq.(3.23) in [14| for complex scalar 
fields] , 

S off = -Trp eff In p cS = ^ [(A/k + 1) ln(M k + 1) - M= In A/" k ] , (78) 

k 

which is also valid for all values of A/k and is much less than the vN entropy in n^. Again, as discussed in Sec. IIV A[ 
the value of the above S e & is the same as the value of the entanglement entropy between the adiabatic particles with 
k and — k, while S e g will decrease as more and more phase information is resolved and the off-diagonal elements 
manifest. 

Unfortunately (|76|) is not well defined if fit = 0, which occurs in the case m 2 /H 2 < 2 when the physical wavelength 
of the mode crosses the horizon. At that very moment, A/"k diverges and the WKB approximation (f76|) fails. In 
this case, after \r]\ = (2 — (m 2 / H 2 )) / k 2 , when the wavelength of the field mode is longer than the size of the horizon 
(superhorizon), the notion of adiabatic particle is no longer viable until the inflation era ends and the universe becomes 
radiation dominated. 



VI. SUMMARY REMARKS 

We conclude with two remarks pointing to the main themes stated in the Introduction, namely, quantum en- 
tanglement depends on partition as well as the choice of physical variables or measurables. They pertain to 1) the 
conditions of separability of quantum states in particle creation processes and 2) the relation of entanglement dynamics 
and entropy generation both measured by the von Neumann entropy. 

A. Conditions of separability of quantum states 

For the model of a free scalar field theory in a dynamical background field or spacetime we see clearly that different 
quantities imply different separabilities of quantum states in different partition and different measurables. Using the 
vN entropy as a common currency for comparisons, we see that 

1. The vN entropy of p = U k + p k = IIk + *k*k with §M for ( k : ~ k ) mode 

pair, or of the Wigner functions (|32p 

or (|34]l. vanishes. This means that p is a pure state, and the quantum field here is a completely isolated system even 
in a classical dynamical background field or spacetime. The vanishing vN entropy of (|45p and f)47[) is similar. 

2. The vN entropy of Tr p ^± k {p} for some specific k with p = nk +v ^k*k' 62)' or (EK i s zero. This means that 
each (k, — k) mode pair is separable from other mode pairs, and of course, each \l/k is pure. 

3. The vN entropy of the reduced Wigner function obtained by tracing out, say, the k 3 < components of (|32p . 
is zero. This means that the quantum state of 5 k 5 k is separable from i?_ k i?I k and all other field modes, while no 

particle with ±k associated with the particle counter or the number operator defined by & k (f?) and 6^(77) is created. 
The zero vN entropy of the reduced Wigner function from (|45p is similar. 

4. The vN entropy (f39| of (|37|) is nonzero, meaning that the particles with k associated with the particle counter 
defined by & k (f?o) and & k (?7o) are entangled with their — k partners. The nonzero vN entropy of (|4"5f has a similar 
meaning. 

5. The vN entropy of the exact density matrix of the vacuum p is zero, but the vN entropy of the effective density 
matrix in Fock representation with off-diagonal elements averaged out is not. The latter indicates that the phase 
information has been coarse grained. 

6. The value of the vN entropy of the effective density matrix (|6ip with all off-diagonal elements averaged out [see 
statements above (|59)) ] is the same as the value of the vN entropy of the reduced density matrix (l6Cfl) after tracing 
out the particles with — k. This suggests that the tracing-out process in obtaining the RDM represents the strongest 
coarse graining. 

7. In the adiabatic number basis, the vN entropy (175)) describing the entanglement between the adiabatic particles 
with k and their — k partners has a different value from (|39p . due to a different choice of physical measurables. 



14 



B. Entanglement entropy and entropy in statistical mechanics 



As we saw in the above the von Neumann entropy has also been used as a measure of the entropy generation in 
particle creation processes. In [l4| KME argued that the effective density matrix p cfi would appear as a mixed state 
when the off-diagonal elements oscillate too rapidly to be resolved. We see close similarity between the vN entropy 
of a bosonic field used in this context of nonequilibrium mechanics and that measuring the quantum entanglement 
between the (k, — k) particles in a single mode pair. Indeed, from our result we see that the correlation between 
each particle pair seems to be equal, so that the entanglement entropy between k and — k particles seems to be 
counting the number of the degrees of freedom that are integrated out. Although the entropy of the former has the 
same value as the latter in the cases considered in this paper, the differences between these two forms of entropy are 
perhaps more revealing, especially when viewed from their respective theoretical structures. 

Let us compare the difference between the vN entropy S = — Tip log p [31] of a closed quantum system with 
the Boltzmann's entropy in a microcanonical ensemble 5* = — fcslogfi where fi is the number of accessible states. 
Note that both describes an isolated quantum system. The enumeration of accessible states is independent of the 
representation and can contain both entangled states and separable states. When one assumes that = Tip is 
given by the probability Trp of finding an isolated system in a particular quantum state, one has already ignored 
all physical information contained in the off-diagonal components of the density matrix such as quantum phase of 
quantum states. This is accomplished under the random phase approximation whereby one can use the concepts of 
probabilities exclusively to describe all statistical mechanical properties of the system. When a system can occupy 
all of its accessible states with equal a priori probability, then the system is in equilibrium. These are, as we know, 
the two fundamental postulates of equilibrium statistical mechanics, namely, (1) equal a priori probability, and (2) 
random phase [32j . (For a depiction of how a quantum system in contact with a thermal bath turn from a quantum 
fluctuation dominated phase to a thermal fluctuation dominated phase at increasing temperature, and under what 
conditions will these two postulates be satisfied in an open quantum system, see [33j.) Finally it is when the particle 
number % > 1 that one can begin to use thermodynamic arguments. 

Thus we see the three stages distinctly: vN and Boltzmann entropy (in its original form) both describe fully isolated 
quantum systems. When one begins to use probability for the description of the system, quantum phase information 
is lost. When one imposes in addition the equal a priori probability assumption one reaches an equilibrium condition. 
Thermodynamic description requires an additional assumption that both the number of particles and the volume of 
the system approach infinity while their ratio is kept a constant. As is known and shown in some of our earlier work 
[jyj the thermodynamic entropy is different from the (equilibrium) statistical mechanical entropy and the quantum 
(nonequilibrium) vN entropy, in increasing order of specificity. 

In contrast, it is clear that the von Neumann entropy of one of the two parties of an isolated, bipartite system is a 
well-defined measure of entanglement [l|| and (|3"5)) is a good entanglement entropy for all for all time. 
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